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Abstract 

We compute the 1-loop chiral corrections to the bag parameters which are needed 
for the discussion of the SUSY — mixing problem in both finite and infinite 
volume. We then show how the bag parameters can be combined among themselves 
and with some auxiliary quantities and thus sensibly reduce the systematic errors due 
to chiral extrapolations as well as those due to finite volume artefacts that are present 
in the results obtained from lattice QCD. We also show that in some cases these 
advantages remain as such even after including the 2-loop chiral corrections. Similar 
discussion is also made for the — > vr electro-weak penguin operators. 



1 Introduction 



We are entering the era of the large scale unquenched numerical simulations of QCD on the 
lattice and so the error on the Bk parameter, currently dominated by the systematic error 
due to quenching [1], is likely to fall below 10% quite soon. This will further improve our 
knowledge on the shape of the CKM unitarity triangle [2], i.e., on the value of the CKM 
phase which is responsible for all the CP- violating phenomena in the Standard Model (SM). 
Since the CP- violation in SM is too small to explain the dynamical generation of the baryon 
asymmetry of the Universe [3] , one is tempted to look for additional sources of CP-violation 
beyond the SM. Supersymmetric (SUSY) extensions of the SM, besides providing an elegant 
solution to the hierarchy problem, also provide new CP-violating phases whose size can 
be constrained by the experimentally measured processes governed by the flavor changing 
neutral currents (FCNC). A convenient way to study those is by using the mass insertion 
approximation [4]. In the basis in which the couplings of quarks and squarks to the neutral 
gauginos are flavor diagonal, the flavor changing insertions arise from the small off-diagonal 
terms in the squark masses, parameterized by dimensionless complex parameters 



y% MB - —2 , (IJ 



where rriq is the diagonal squark mass (averaged as rriq = ^Jfn^lrT^), and {/\^-'^)ab are 
the off-diagonal elements which mix both the left-handed and right-handed squark flavors 
(i, j = 1, 2, 3; A, 5 = L, R). Like in the SM, Sk, the parameter which measures the indirect 
CP-violation in the neutral kaon system, is given by 

but unlike in the SM, the effective Hamiltonian in SUSY [H^^^"^ = J2i Ci{fi)Oi{n)], besides 
the left-left {LL) four-quark operator, also contains the LR, RL and RR ones (see below 
for the specific bases of such operators). The Wilson coefficients, Ci{fi), can be computed 
perturbatively in any specific low energy SUSY model. The matrix elements, instead, must 
be computed non-perturbatively. To date the most suitable tool to do such a computation is 
by means of lattice QCD. The matrix elements of the non-SM operators are enhanced with 
respect to the SM one by a large factor, m\/{ms + m^)^ 25, and therefore their more 
accurate determination is mandatory. So far these matrix elements have been computed 
in quenched lattice QCD by considering kaons consisting of degenerate quarks [5]. The 
corresponding results are then used to discuss the constraints on (5^2)^5 [6], along the lines 
proposed in ref. [7]. In addition to unquenching, the lattice estimates should be improved by 
fixing one of the kaons' valence quarks to the physical strange quark mass (accessible in the 
lattice simulations) and pushing the other one as close to the chiral limit as possible. Due to 
limited computing resources, however, two problems appear: (1) one cannot work with the 
light quark as light as the physical (i-quark and therefore a chiral extrapolation will always 
be needed; (2) the finite volume effects become more pronounced as the light quark mass 
is decreased. In such a situation and in order to get the phenomenologically interesting 
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results from lattice QCD on the matrix elements relevant to the SUSY — mixing 
problem, one should find a way to reduce uncertainties related to these two problems. In this 
paper these issues are addressed and further considered by using chiral perturbation theory 
(ChPT). We compute the chiral logarithmic corrections to the so called "bag" -parameters 
and discuss the possible strategies that would allow one to minimize their impact onto the 
chiral extrapolations of the lattice results, as well as to minimize the systematic errors due 
to finite volume. In Sec. 2 we recall the frequently employed bases of AS = 2 operators and 
define the bag-parameters, to which we compute the chiral corrections in Sec. 3. In Sec. 4 
we discuss the combinations of bag parameters which are (completely or partially) free of 
chiral logarithms, and we briefiy comment (Sec. 5) on the application of a similar strategy 
to the K —>■ 71 matrix elements that are relevant to the A J = 3/2 amplitude in K irir 
decay. In Sec. 6 we discuss the impact of the 2-loop chiral corrections to one class of the 
combinations of bag parameters, and we briefiy conclude in Sec. 7. 



2 Bases of AS = 2 operators and parameters 

The SUSY contributions to the K^ — K^ mixing amplitude are usually discussed in the so 
called SUSY basis of AS* = 2 operators [7]: 



Oi = 






02 = 


s\l - 




03 = 


s\l - 




04 = 


s\l - 


75)rf*s^■(l+75)c^^ 


05 = 


s\l - 


75)rf^s^'(l+75)c?% 



where the superscripts stand for the color indices. Other bases are also employed, of which 
the Dirac basis is probably the most popular among the lattice QCD practitioners, 

gl = s*7^(l-75)rf*s^7M(l-75)rf^ 

Qz = s\l + ^,)d'7s^{l-^,)d^ , (4) 

g4 = s*(l-75)^^'s^■(l-75)c^^ 
= sV^^rf* Pa^yd^ , 

where we use the definition in which a^y = [7^,7i,]/2. Although the operators in the above 
bases are written with both parity even and parity odd parts, only the parity even ones 
survive in the kaon matrix elements. The latter are usually parameterized in terms of 
bag-parameters -B1-5, namely. 
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(if"|04M|ii-»> = 2 ( , )'r4/i.i?,M. 

3 \ms{u) + md[v)J 

where 5^ > 0. z/ in the above equation indicates the renormahzation scale of the log- 
arithmically divergent operators, Oj, and the scale at which the separation between the 
long-distance (matrix elements) and short- distance (Wilson coefficients) physics is made. 
To make contact between the matrix elements of the operators in (3) and those in (4), one 
applies the Fierz identity on the Dirac structures (FD), which leaves the physical amplitude 
invariant, and then reorder the color reversed indices in operators O3 and O5. For example, 

— {Sadls'A [1® 1+ 75 ® 75 - ® ^f^uU,^^) = \ ((Qs) - (Q4)) . (6) 
Similarly, (O5) = —(1/2) ((52)- Summarizing, 

(Qi) = (Oi) , {Q2) = -2(05) , 



(g5) = (02)+2(03). (7) 

In eqs. (6,7) and in the rest of this paper the i^-dependence is implicit. As a side remark, 
we note that these formulae are strictly true only in the renormahzation schemes in which 
the Dirac Fierz identity is not violated, such as the (Landau) RI/MOM scheme [8]. ^ In the 
following we will suppose that the subtraction of ultraviolet divergences is made in such a 
renormahzation scheme and restrict our attention to the low energy behavior of the above 
matrix elements. 



3 Chiral logarithmic corrections 

As mentioned in introduction, we will use ChPT to discuss the low energy behavior of 
AS = 2 operators relevant to the SUSY K^-W mixing amplitude. Even before considering 
the chiral representation of eq. (3) or eq. (4), it is clear that the chiral behavior of the 
matrix elements of O2 and O3 will be the same since these two operators differ only in the 

"'^See ref. [9] for a formulation of the MS scheme in which FD is preserved at the next-to-leading order 
(NLO) in perturbation theory. 
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color indices. In other words these two operators differ by a gluon exchange, which is a 
local effect, that cannot inffuence the long distance behavior described by ChPT. Therefore 
the chiral logarithms in B2 and B3 will be the same although their respective low energy 
constants (LEC's) are different. The same argument applies to B4 and B^. This "color 
blindness" is evident when working out the chiral representation of the operators (4). To 
that end, we will use the lagrangian and notation specified in our previous paper [10], and 
account for the following properties: 

- Under SU(3)l ® SU(3)r the field S transforms as S — RT,L^; 

- The lowest order Lorentz scalars, transforming as (3,3) and (3,3), are E and E''', 
respectively; 



- The lowest order Lorentz vectors, transforming as (8,1) and (1,8), are zE'''(9^E and 
iTid^Ti^ respectively. 

We are now able to write the bosonised versions of eq. (4), namely. 



Qi 
Q2 
Q3 
Q4 
Q5 



ds 



(8) 
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'^ds^ds 



^ds^ds 



"^ds^ds 



^ds^ds 



where we introduced the new set of bag parameters, 6j, with the signs chosen as to make 
all 6j's positive. After sandwiching the above operators by {K^\ and \K^) and after eval- 
uating the matrix elements at leading order, we can relate 6j's to the chiral limit of the 
5i-parameters: 
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(9) 



We now proceed, like in ref. [10], by following the standard routine to compute the chiral 
logarithmic corrections to -Bj. In computation of the tadpole chiral loop integrals we use 
the naive dimensional regularization and the MS renormalization scheme. Our results are: 
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(10) 
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where dots stand for analytic and higher order terms in ChPT, and /i is the renormalization 
scale. 



4 Log-safe combinations 

The phenomenological applications of the predictions based on ChPT at NLO are usually 
plagued by the poor knowledge of the size of low energy constants [11]. A better pre- 
dictability is then expected for the combinations of physical quantities in which the low 
energy constants cancel (partially or completely). Contrary to that situation, when com- 
puting the physical quantities from the QCD simulations on the lattice one works with light 
quark masses larger than the physical up- or down-quark {rriu/d), thus allowing one to probe 
the analytic dependence on the quark masses, while missing (again, partially or completely) 
the chiral logarithmic behavior that is expected to take over as the light quark becomes 
closer to physical rriu/d- Since the point at which the chiral logarithms, with coefficients 
predicted by 1-loop ChPT, are to be included in extrapolations of the lattice results is not 
known, their inclusion in the chiral extrapolations induce large systematic uncertainties. To 
avoid such uncertainties one should aim at combining the physical quantities in which the 
chiral log corrections cancel. As we will show, such log-safe combinations also help reducing 
the finite volume artifacts that are becoming ever more important as the light quark gets 
closer to the chiral limit. 

We now construct the golden (silver) combinations in which the chiral logarithms com- 
pletely (partially) cancel. The criterion for creating the silver combinations will be the 
cancellation of the pion loop sum or integrals because they make the strongest deviation 
in the chiral extrapolations of the results obtained directly from lattice QCD, and because 
they generally represent the most important source of the finite volume artifacts. From the 
discussion in Sec. 3 and eq. (10), one immediately identifies the following two golden ratios: 

As for the silver combinations the simplest ones that can be deduced after inspecting eq. (10) 
are 

Rtj = B,xBj (^ = 2,3; J =4, 5). (12) 

Alternatively one can use some auxiliary quantities, preferably those that are easily cal- 
culable on the lattice, and combine them with bag parameters (10) to cancel the pionic 
logarithms. Useful quantities are the decay constants and their combinations [12]: 
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(13) 
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The silver log-safe combinations in which the decay constants are combined with i?-parameters 
are: 



f2 

Jk n 
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72"^4,5 

Jk 



(14) 



If, as widely expected, the i?i 
log-safe quantities are also 



Rp. 



Bk parameter is accurately determined first, other silver 

-^2,3 



5i 



Rj = -84^5 X Bi 



(15) 



Summarizing, in the golden ratios, R121 all chiral logarithms cancel, whereas in the silver 
combinations, R^^ and Rf^j, the cancellation of the most problematic part (from the lattice 
practitioner's point of view) is achieved. 



4.1 Finite volume effects 



In the finite box of side L, instead of integrals one deals with the sums over discretised 
momenta q = (27r/L) x n [n G X^]. The difference between sums and integrals is the 
infrared (/i independent) effect that can be expressed in terms of the function "^s" whose 
properties we discussed in ref. [10]. In other words, for large physical volumes, one can 
deduce the finite volume effects by comparing the expressions for a given physical quantity 
derived in ChPT in finite and infinite volumes [13]. Like in ref. [10], we define 

Bp 



ABi 



B?" 



B^ 



(16) 



and obtain: 
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AE4,5 
-84,5 



2/2 

1 

2f 



2/2 



m 



K + '^i, (T ^ ^ 
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771^^3 {L,mK] 



7m^ 



m 



K 



2 O ^ 



:i7) 



-a (L, m^) + 2m^^3 (L, tuk) + -6 {L, m^) 
z ^ 2 ^ 



It is important to notice that the coefficients multiplying the function ^1 ^{L^mp) are the 
same as those in eq. (10) multiplying mp log(mp//i2). ^ Therefore the combinations of 

^The factor 2 of mismatch between the coefficients multiplying the kaon part in (10) and in (17) is 
canceled by a factor of two in the definition of ^3 (L, tok). 
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physical quantities in which the chiral logarithms cancel not only allow for the safer chiral 
extrapolations of their lattice estimates, but they also provide the cancellation of the finite 
volume effects (at least those predicted by the 1-loop ChPT). For the silver combinations 
there is however a subtlety: although exponentially suppressed [oc exp(— m/^L)], the terms 
proportional to ^3(L,mx) are numerically important because they have larger factors in 
front and because the function ^3 has worse infrared behavior than ^1. Therefore it is 
also important to be careful in dealing with terms corresponding to the Icaon loops. Of all 
the silver combinations discussed in this section, only R^j [see eq. (12)] receive large finite 
volume corrections while all the other silver combinations do not suffer from this problem. 
The finite volume effects for RI_y, together with those for -Bi_5, are plotted in fig. 1. We 
first observe the usual behavior, namely that the finite volume effects become larger as one 
is getting closer to the physical point, i.e., r = 0.04 [14]. Secondly, we see that the finite 
volume effects on the silver combinations, RI_y, are clearly reduced when confronted to 
those that plague i?i_5-parameters. Finally, we stress again that the leading finite volume 
effects on the golden ratios [eq. (11)] are totally absent. 

5 K ^ 7T case 

The discussion of the previous section can be easily extended to the K —>■ n case which 
is often considered on the lattice when computing the A/ = 3/2 amplitude in K tttt 
decay. A golden log-free combination can be easily constructed from the chiral logarithmic 
corrections to the electro-weak penguin operators calculated in, for example, ref. [15]. 
To be more specific we will concentrate on the following three operators: 
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where 7^^ 



R/L 



7^(1 ± 75). The chiral representation of these operators reads 




(19) 



The relevant matrix elements are parameterized as 
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Figure 1: Finite volume effects for the Bi-^-paramctcrs, computed by using eq. (17), are denoted by 
the thick sohd curves, whereas those corresponding to the silver combinations i?f _7 are represented by the 
dashed and dotted thick lines. The plots are made for L = 2fm. Notice that r = m^/d/'nT's is varied by 
keeping mg fixed to its physical value. To better appreciate the benefit of using R\_-j we also show the line 
corresponding to zero finite volume effect (thin dashed line). 



where the bag parameters are the K ^ n versions of the — ones, discussed 

in the previous section. They both have the same tree level values but differ in the chiral 
logarithmic corrections, which we computed as well and obtained 



tree 



log + 2 tjik 



4(47r/)2 \y rriK — rn-j^ ^ " /^^ "^k — 



log 



Kit 



b: 



tree 



1 + 



log — 

rriK — niT, fi- 



-3m^log-^J + 

rriK - 5m^ 
h 2 m 



TTlK 



,2 



2 i_ 

K 



log- 



7 2, "^r) \ 

-3m,log-)+... 



(21) 



^5 



Dtrcc 
^5 



4(47r/)^ 



m: 



m^log^ + 2 



-m 



2 



+ 



log 



Therefore, the chiral extrapolation of the ratio of matrix elements of the electro-weak pen- 
guin operators, i.e.. 



d9 K-k 

JXry 



(22) 



that can be computed on the lattice, is free from uncertainties induced by the inclusion of 
the chiral logarithms in chiral extrapolations. In addition, the leading finite volume effects 
in the ratio -Rf ^^'^ completely cancel. The absolute values for -Bf/5, on the other hand, can 
be extracted from the following silver relations: 



s Ktt 



R. 



s Kit 

■4,5 



f tdKtt 



(23) 



We believe the safe extrapolation of the lattice results for may be useful in disentan- 

gling the current discrepancies among various analytic [16] and lattice approaches [17] used 
so far to estimate {7r^\Qg^'^\K~^) . Even more so after noticing that almost all approaches 
agree in the value for B^'^, while they differ quite a lot in that for B^'". Finally, a reader 



interested in the discussion of the potential SUSY enhancement of the AI 
is referred to ref. [18]. 



3/2 amplitude 



6 Golden ratios and next-to— next-to leading (NNLO) 
chiral corrections 

We showed that the golden ratios were completely protected against the presence of the 
chiral logarithms and thus also against the finite volume corrections. One may wonder if 
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that feature survives after accounting for higher order corrections in the chiral expansion. 
That question was recently rised in ref. [19] where it has been argued that the pion mass 
might receive sizable finite volume NNLO chiral corrections. That issue will obviously be 
settled only after the full 2-loop chiral corrections in finite volume are computed. Here we 
discuss such corrections to our golden ratios -Rf 2, and argue that in this case the finite 
volume corrections arising beyond 1-loop ChPT are indeed negligible. Golden ratios are 
generically defined as 

^ {out\OA\tn) 

{out\OB\in) ' ^ ' 

where Oa and Ob are two different operators with the same chiral representation. The 
NNLO chiral expansion to this ratio can be schematically written as 

^ 5'a[1 + {logs + ca) + {dlogs + clogs a + Ca) + • • •] ^25) 
qbII + {logs + cb) + {dlogs + clogs b + Cb) + • • •] 



9b 



|l + (c^ - Cb) + [{dogs A - clogsB) + {cb - ca){cb + logs) + {Ca - Cb)] + •••}, 



where gA/9B is the tree level value of the golden ratio, ca,b and Ca,b are the low energy 
constants arising at NLO and NNLO respectively, clogs a,b are the 1-loop part to the NNLO 
correction. Based on the same argument used in sec. 3, the single and double chiral log- 
arithms, logs and dlogs, completely cancel in the ratio and thus the computation of the 
NNLO chiral corrections comprises the 1-loop diagrams only. That feature clearly holds at 
any other higher order, i.e. one needs to compute the diagrams that are one loop less with 
respect to the ones indicated by the canonical ChPT counting. 

After including the NNLO corrections, the finite volume effects to the golden ratio can 
be written in the following form: 

= = AclogsA - AclogsB + {cb - CA)Alogs , (26) 

in an obvious notation. We see that the finite volume effects arise only from the 1-loop 
diagrams at the weak vertices, and therefore the result can be expressed in terms of the 
^^-function, as before. 

To exemplify the above discussion we now use the set of counterterm lagrangians, written 
explicitely in ref. [20], to obtain the counterterm contributions to the NLO expression for 
Those are then used to calculate the finite volume corrections to i?f at NNLO. We 
obtain ^ 

,g = -j^{l2ml [5i + 3^2 + 2 (-353 + 8(54 + 8^5 + 856 -457)] ^{L,mK) 



ARl 



Rl 12/2 



■^To easily identify the counterterm coefficients from ref. [20], our (5's and their c's are related as 

94 95 

with .94 = 23^/^31'"", and 55 = {2/ 3)3^/^31""". 



10 



+3 [8m%62-m^^ {6i + 62 + i (363-56^-565-66^ + 667))] ^{L,m^) 
+ [32m% (^4 + ^5 + Sq) - 4ml {S6^ + 36^ + 2 6^) 

-m\ (5i + 17 ^2 + 20 53 + 24 ^7)] ^ (L, m^) } . (27) 

Notice that the double pole contributions (oc ^3) coming from NLO^ and from NNLO terms, 
cancel against each other. The size of the finite volume corrections evidently depends on 
the values for the low energy constant 6i. By using either the naive dimensional analysis or 
the values extracted from the recent quenched lattice data ^, we get that for L > 1.5 fm 
and r G [0.04, 1), the finite volume effects, AR2/R2, are within a few per mil, thus totally 
negligible. 

7 Summary 

In this work we computed the 1-loop chiral corrections to the bag parameters that are 
relevant to the K'^ — mixing amplitude in the SUSY extensions of the SM. Two out 
of five independent AS* = 2 operators have reversed color indices so that their chiral log 
corrections are identical to those in which the color indices are not reversed. This "color 
blindness" , which also emerges from the explicit calculation in the chiral representation of 
those AS* = 2 operators, is actually expected since the change of the color structure of 
the operators is a local effect which cannot appear in chiral loops. Rather it is encoded 
in the low energy constants. After comparing the 1-loop chiral expressions for all the bag 
parameters obtained in finite and infinite volume, we also provide the formulae for the finite 
volume effects which are useful for the assessment of the associated systematic uncertainty 
in the lattice calculations. As in our previous paper [10], we see that the finite volume 
effects become more pronounced as the light valence quark of the neutral kaon gets closer 
to the chiral limit. We find that the combinations of bag parameters in which the chiral 
logarithmic dependence cancels partly or totally may be particularly beneficiary for the 
lattice calculations because they allow one to sensibly reduce two important sources of 
systematic uncertainties: (i) the errors induced by the chiral extrapolations, and (ii) the 
errors due to the finiteness of the lattice box. We construct the explicit combinations of 
bag parameters that are log-safe, i.e., in which the chiral logarithms cancel either totally 
(golden combinations, -R12) or partially (silver combinations, Rl_i)- In addition, we show 
that inclusion of the NNLO corrections does not spoil the advanteges of computing the 
golden ratio on the lattice. A similar discussion is also made for the case of the K ^ tt 
matrix elements of the electro-weak penguin operators. 



"'We thank Mauro Papinutto for sending us the lattice estimates for (J^'s prior to their pubhcation. 
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